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Abstract 

This paper tells the story how a MAPLE-assisted quest for an interest- 
ing undergraduate problem in trigonometric series led some "amateurs" 
to the discovery that the one-parameter family of deterministic trigono- 
metric series S p t \-t ^2 n& ^sm(n^' p t), p > 1, exhibits both order and 
apparent chaos, and how this has prompted some professionals to of- 
fer their expert insights. As to order, an elementary (undergraduate) 
proof is given that S p (t) = a p sign(i)|i| 1 /P + 0(|<| 1 /(p+ 1 )) V t £ R, with 
explicitly computed constant a p . As to chaos, the seemingly erratic fluc- 
tuations about this overall trend are discussed. Experts' commentaries 
are reproduced as to why the fluctuations of S p (t) — a p sign(i)|t| 1//p are 
presumably not Gaussian. Inspired by a central limit type theorem of 
Marc Kac, a well-motivated conjecture is formulated to the effect that 
the fluctuations of the [i 1// ^ p+1 ^]-th partial sum of S p (t), when properly 
scaled, do converge in distribution to a standard Gaussian when t — > oo, 
though — provided that p is chosen so that the frequencies \n p }nsN 
are rationally linear independent; no conjecture has been forthcoming 
for rationally dependent {n _p } ne psj. Moreover, following other experts' 
tip-offs, the interesting relationship of the asymptotics of S p (t) to prop- 
erties of the Riemann £ function is exhibited using the Mellin transform. 
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1 Introduction 



Back in the 1990s when I was one of Jerry Percus' postdocs, I learned that 
Jerry's curiosity often let him explore unorthodox scientific ideas, just to see 
where they would lead to. In this vein, I take the invitation to celebrate the 
seminal contributions to statistical physics by three of its living legends: Ben 
and Jerry, and Michael, as a wonderful opportunity for me to follow Jerry's 
example and to take the three honorees (and the reader) on a curious trip 
into the realm of deterministic chaos without pretending that I am motivated 
by a physics problem I am not! Neither do I claim any mathematical 
sophistication! It is just an amusing story to tell, involving several actors, 
interesting mathematics, a few rigorous results, and some conjectures. 
The object of study is the one-parameter family of sine series 

s P (t) = E„ eN sin (^); ®p>i, (i) 

which converges absolutely for t G R\ it's not in [Zyg02|. Since t -1 5 p (t)^K (p), 
which is Riemann's Zeta function (art. VII in |Rie76j ) . the study of p t- >■ S p (t) 
for fixed t, when analytically extended to p G C\{1}, might be of interest 
to analytic number theorists. However, I don't know whether this produces 
anything not already known about (. Indeed, some relationships between 
S p (t) and C(s) which go beyond the obvious one just exhibited were pointed 
out to me by Norm Frankel and, independently, Steve Miller, in response to 
my SMM 106 talk. Prompted by their insights I added section 4.2. 

For the most part our attention will be on the t-dependence of S p (t) for 
real p > 1. Since S p (— t) = — 5 p (t), it suffices to discuss 5 p (t) for positive t. 

Since the sine function with the shortest wavelength contained in S p (t) is 
sin(t), to which sine functions with ever longer wavelengths are being added, it 
is to be expected that S p (t) is neither periodic nor quasi-periodic. Interestingly 
enough, the deterministic map t i— > S p {t) exhibits apparent chaos on small 
scales, yet order on large ones. For example, here are two plots of S2 {t)'- 
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Fig.l. The 5,000-th partial sum of S 2 (t) versus t for < t < 500. 
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The second one is over a 100 times larger interval of t values: 
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Fig.2. The 200,000-th partial sum of J 2 (t) versus t for < t < 50, 000. 

We see that relative to the range of values taken by ^(i), the seemingly erratic 
oscillations around their local mean appear to decrease with increasing S2 value 
range, and the graph appears to converge onto a rightward opening parabola, 
an increase of the t domain by a factor 100 producing an increase of the J>2 
value range by a factor of 10; i.e., a square root type behavior. Qualitatively 
similar p-th root type trends of S p {t) can be observed for other values of p > 1. 

One of the rigorous results to be proved in this paper, with elementary 
means, is that S p (t) = a p t l / p + 0(t 1/(p+1) ) V t > 0, with explicitly determined 
a p for all p > 1. This was obtained in partial collaboration with Jared Speck. 

For small t, numerical evidence is given that an 0(t 1 /( p+1 )) bound on the 
deviations from the overall trend is optimal, while it becomes lousy for large 
times. As pointed out by one of the three referees]]] improved bounds on the 
deviations from the trend for large t can be obtained if the Riemann hypothesis 
is assumed. I summarize their comments in the added section 4.1.1. 

More difficult than the determination of the trend function, but also more 
interesting, is the analysis of the deterministic, yet apparently chaotic fluctu- 
ations about the overall trend. A discussion of Kac's central limit theorem for 
sine series with rationally independent frequencies will lead us to the conjecture 
that the fluctuations of the |t 1 '( p+1 )]-th partial sum of S p (t), when properly 
scaled, do converge in distribution to a standard Gaussian when t — > 00 - 
provided that p is chosen so that the frequencies {n~ p } ne ?q are rationally lin- 
ear independent]^] no conjecture has been forthcoming for rationally dependent 

1 Any similarity with the number of honorees is unintended and purely coincidental. 

2 It is clear that p must be chosen irrational. However, as noted by one of the referees, 
p Q is not sufficient to obtain rationally linear independent frequencies of the form rT v : 
namely, the set {n _p } nS N will be rationally linear dependent whenever p = with integers 
a > b > 1, and this formula produces rational as well as irrational p. 
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{n~ p } n m- The stronger conjecture that S p (t) — a p t 1 ^ p exhibits Gaussian fluc- 
tuations, entertained by me at the time of SMM 106, is presumably wrong, 
as pointed out to me by two of the expert referees; see section 4.1.1. Perhaps 
the discussion will prompt some interested reader to work out the definitive 
answer using the professionals' tools. 

Before we now plunge into the rigorous analysis of the functions t r-> S p (t), 

1 owe the reader an answer to the burning question: How come I got to dabble 
in the math of these sine series? After all, this is not my field of expertise! 
The answer is: A question by my colleague Steve Greenfield about the graph 
of Jhit) for < t < 120 originally got me started, and the rest was curiosity 
about the behavior of ^(t) for later t, and some fascination with what I found. 
So I begin with the simpler (but not so simple) behavior of S p (t) at early times. 

2 The early time behavior of S p (t) 
2.1 "Can you explain the tilt?" 

On April 10, 2007, Herr Dr. Prof, (emeritus) Stephen Greenfield^ sent me the 
following email: 

"The attached picture is a graph of the 100th partial sum of the infinite 
series whose nth term is sin(x/n 2 ). You are a clever person. Why does 
the graph have the "tilt" that it does?" 




— i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i 

O 25 50 75 lOO 



Fig. 3. (Color online) Steve Greenfield's Maple plot of S2 (%)■ 

Nice question; but before I would drop everything I was doing at the time to 
rise to the challenge, I wanted to know why he was looking at that trigono- 
metric series. So I went to Steve's office three doors down the hallway to ask 

3 Steve likes to make fun of a German convention by calling me "Herr Dr. Prof. Kiessling;" 
so I assume it's only fair when I reciprocate. 
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him what this was all about. As it turned out, while trying to invent some 
interesting unorthodox calculus problem for his honors undergraduate Maple 
workshop, one you won't easily find solved in a solutions manual, he had played 
with some unconventional trigonometric series, and this one exhibited some 
curious behavior: Why does the oscillating graph show some overall upward 
trend, instead of oscillating about zero, like more conventional sine series? Is 
there an explanation which a good undergraduate student could understand? 

An elementary, positive lower bound which tilts upward over the full do- 
main displayed in Steve Greenfield's picture was soon found. By itself this 
bound does not suffice to explain the overall shape of the graph of S2 (t) but 
at least it explained why the graph of S2(t) wasn't oscillating about zero. More 
importantly, however, it would open the flood gates and let curiosity as to the 
behavior of ^(t) take hold of me, and others! This bound is reproduced below. 

Recall that for t > we have sin(t/n 2 ) > t/n 2 — t 3 /6n 6 . For t/n 2 < 1, this 

lower bound of sin(t/n 2 ) is off by 16% at worst. It can be used in the series 

defining ife(t) whenever n > \\/t], where |Y] is the smallest integer not less 

than the real number r. Thus, writing 

Wt\ 00 
S 2 (t) = £ sm(n- 2 t) + £ sin(n- 2 t), (2) 
n=l n=|y*l+l 

we estimate the second sum from below by 

£ sin(n- 2 t)>( £ n- 2 )t-\( £ n^)t\ (3) 

n=\y/t\+l ^ n =\y/t\ + l 7 V n=|yt]+1 7 

Furthermore, by the familiar Riemann sum approximations, we estimate 

oo r°° 

E £>/ ^=T7fei' ^ 



E < / ^^-|t4t5> (5) 
and so we find 

oo 

E sm(n-H) > jj^ - ± oJ ^ . (6) 

R.h.s.([6]) is a piecewise cubic lower bound to l.h.s.((6]), and therefore much 
easier to discuss than l.h.s.(|6]). In particular, it is easily seen to be positive 
and overall increasing roug hly like y/i; see Fig. 5 below. On the other hand, the 
first sum at r.h.s.© has just \y/t \ < 11 terms for the t (viz. x) interval shown 
in Fig. 3. These were few enough to show by direct inspection that it did not 
have enough negative terms to overpower r.h.s.([6]). In an undergraduate class 



I am resorting to my choice of variable "£" rather than Greenfield's "x" because the 
overall thrust of my paper is to think of 1 1— > S p (t) as a deterministic process in time. 



5 



one would simply have to allude to the fact that a finite sum with less than a 
dozen terms is manageable and not go into details, though^ 

Although I didn't plot it at the time, a Maple plot of the first sum at 
r.h.s.([2]) shows that it itself is non-negative and tilted upward for < t < 120: 




Fig. 4. The first sum at r.h.s.JJ}. 

So, curiously enough, Steve Greenfield's question applies verbatim to Fig. 4! 
Yet, rather than trying to explain the overall upward tilt in Fig.4, one may 
want to try to prove merely that the first sum at r.h.s.Q is non-negative for 
< t < 120. I haven't tried it, but the upshot of any such proof is: r.h.s. fifty 
is an elementary lower bound to for < t < 120. This bound is actually 
quite decent; see Fig. 5 below. 




Fig.5. The 2,000th partial sum of S2(t) together with r.h.s. JBJ. 



5 A suitable "undergraduate bound" on 52 (t) is supplied in the appendix, however. 
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2.2 Do all series S p (t) have graphs like that of S2 (t)? 

We pause briefly to inspect the early time behavior (up to t — 120) of some 
sine series with other parameter values p > 1. Here are a few examples. 
The first figure shows the graph of S p (t) for p = a/2: 

so- / 
50 
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t 

Fig.6. The 300,000-th partial sum of S^{t). 

That graph looks comparable to that of 5 2 (t), only that the overall tilt is 
steeper, roughly by a factor six. The amplitudes of the oscillations in the 
graph of 5^/2(2) appear smaller than in the graph of S2(t), but appearances 
are misleading, for the overall range of values is about six times as large. 
In absolute terms the local fluctuations actually have increased, from a local 
amplitude of 1-2 in the graph of ^(i) to about 2-3 in the graph of S^. By 
the way, by "amplitude" I mean half the difference between a local maximum 
and its ensuing local minimum in the graph. 

Next we see the graph of S p (t) for p = \fl . It is plotted separately because 
it would show merely as a "bottom dweller" if incorporated in Fig.6. 
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Fig. 7. The 2,000-th partial sum of S^(t). 
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Also this graph looks comparable to that of S2 (t). Now the overall tilt is less 
steep, roughly by a factor one half. The amplitudes of the oscillations in the 
graph of Jyf (t) appear larger than in the graph of S2 (t), but again appearances 
are misleading, for the overall range of S^m is about half as large. In absolute 
terms the local fluctuations have decreased, from a local amplitude of 1-2 in 
the graph of S2 (t) to something closer to 1 in the graph of 5^- 

Note that the oscillations about the upward trend show GO/tt ~ 20 minima 
in all three figures, corresponding to the shortest wavelength involved. 

Lest the reader now thinks that, except for the magnitude of their tilts, 
the graphs of S p (t) look roughly alike for all values of p > 1, we note that 

lim S p (t) = sin(t), ViGf. (7) 

Therefore, eventually the graph of S p {t) will look essentially like that of sin(t) 
over the whole t interval [0, 120]. (I spare the reader the graph of sin(i).) 

The discussion in the previous paragraph also makes it plain that the strict 
positivity of all displayed S p (t) graphs for t > is due to a too small sample of 
p values near p = 2. Eventually when p is large enough, the graph of S p {t) will 
cross the t-axis. By comparing Fig. 3 with Fig. 7 it should come at no surprise 
that p = 3 is large enough; however, I didn't attempt to determine the critical 
p- value at which the first positive solution to S p (t) = occurs, nor do I know 
whether this would be interesting to know. 

We now turn to the more interesting problem of the overall shape of S p {t). 

3 The overall shape of S p (t) 

My tending to Greenfield's question about the "tilt" of S2 (t) had produced the 
lower estimate to l.h.s.()6]) given by r.h.s.()6]); this estimate is bounded below by 
CVi, and x (29/30) y/i for t — > 00. Moreover, in a similar fashion one obtains 
an upper bound l.h.s.(|6])< C'y/t which is asymptotic to y/t for t — > 00. These 
findings implied that l.h.s.(|6])= C'y/i+small corrections. Furthermore, the 
first sum at r.h.s.()2]) was bounded absolutely by C"'\/t and otherwise should 
be responsible for all the fluctuations visible in the plot. So when I presented 
Steve with my lower bound to S2 (t), I also told him about my conjecture that 
S2 (t) = a 2 V / ^+fluctuations for some constant a 2 . 

The conjecture surprised Steve, for Fig. 3 had suggested to him that the 
graph of S2 (t) will continue to grow on average at roughly the same rate as the 
overall tilt visible in Fig. 3. To be fair, there isn't much of an overall concave 
bent of the graph of S2(f) to be seen in Fig.3. Using Maple, a graph of S2(t) 
similar to the one shown in Fig.l was now produced, and compared with y/t. 
It confirmed the a 2 \/t trend; however, a 2 had to be somewhat bigger than 1. 
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3.1 The pursuit of a 2 



Enter Jared Speck, who at the time worked on his Ph.D. thesis research in 
relativity, advised jointly by me and my colleague Shadi Tahvildar-Zadeh, and 
who even may have been Greenfield's TA at the time. When I told him about 
Greenfield's Szit) and my conjecture about its y't-like trend, he didn't exactly 
drop whatever he was doing at the time, but the problem didn't let go of 
him either. By the end of April 11 (midnight, that is...) he had produced a 
conjecture as to what the constant a 2 could be! Jared noted that by boldly 
replacing the sum over n G N with an integral over "c?n" from 1 to oo, followed 
by the variable substitution /z 2 = t/n 2 , one obtains a factor \ft that can be 
pulled in front of the d\i integral, and letting t — > oo in the upper limit of that 
d\x integral, one obtains a candidate for a 2 , namely 



April 12 was spent pondering Jared's bold proposal. 

On the one hand, there was no reason to expect thafj sin(V 2 t)dv was 
an accurate pointwise approximation to 52(f) as t — > 00 because sm(n~ 2 t) hops 
around erratically in the interval [—1, 1] when n h-> n + 1 for small n (and there 
are more and more "small" n as t becomes large), so that one could not allude 
to a Riemann sum approximation. On the other hand, perhaps we could show 
that the difference between -^5 2 (^) and 4= sin(i/~ 2 t)aV would tend to zero, 
so that his conjecture would be correct asymptotically. 

The first impulse was to resort to the splitting (j2J) of the series defining 
Si{t). The already collected evidence that the second term at r.h.s.(T5]) x C"y/t, 
with C" < 1, suggested that all that needed to be done was to prove that the 
first sum at r.h.s.© made another, though smaller, cx sjt contribution, aside 
from yielding a subdominant erratic behavior. But there was an obstacle. 
Using the upper and lower estimates — 1 < sin(£) < 1 produces upper and 
lower bounds ±y/t on the first sum at r.h.s.(T5]) which, while compatible with 
the required oc y/t contribution, aren't good enough. There must be many 
near cancellations in that sum, but a term-by-term discussion, though feasible 
for small t, was of course out of the question for larger t. 

Later that evening I realized that the key was indeed to split the sum 
of S2{t) into two parts, but not as done in the lower estimate given in the 
previous section — instead, one had to split at some n oc |~t 1//3 ] rather than at 
n oc [t 1 / 2 ]. More precisely, with r chosen <n/2, when t is large enough then 
for n > [(2t/r) 1 / 3 ] any two consecutive arguments t/n 2 and t/(n + l) 2 of the 
sine functions would come to lie within a quarter period of sine; furthermore, 

6 To avoid unnecessary confusion, I switch to v rather then "n" to denote the continuous 
integration variable, and leave n to denote the discrete summation variable. 
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with increasing n, for fixed t/r, the consecutive arguments t/n 2 and t/(n+ l) 2 
would be more and more closely spaced. Put differently, for fixed sufficiently 
small r, with increasing t the part of the sum of with n > \(2t/r) 1 ^ 3 ] 
will be an increasingly better Riemann sum approximation of the integral 
Jr(2tA) 1/3 l+i sin(z/- 2 t)<izA Explicitly, if we split 

S 2 (t)= £ sin(n- 2 t)+ £ sin(n- 2 t), (9) 

then, for large t/r, by Riemann sum approximation and substitution /i 2 = t/z/ 2 , 

,tV2 /r(2t/r) i/3 +11 

£ sin(n~ 2 t) « \/t / [i~ 2 smfi 2 d[i. (10) 

n=[(2t/r) 1 /3] + i JO 

For any fixed r the upper limit of the integral at the r.h.s.(!T0|) grows essentially 
oc t 1 / 6 , i.e. it slowly but steadily diverges to oo as t — > oo, and so this integral 
converges to r.h.s.flEl)- Moreover, the first sum in is obviously subdominant. 
Better yet, this erratic term should have plenty of near self-cancellations, and 
if one could show that it vanished on average, then even Jared's replacing of 
by Jj 00 sin(z/~ 2 t)<iz/ could conceivably be vindicated in an average sense. 
In any event, by now I had become convinced that Jared's conjecture for a 2 
was right, and I sent an email to him and Steve detailing my thoughts. 

An hour or so later, but still the same day (almost midnight, again), Jared 
replied with the following email (temporarily we are back to x instead of t)0 

" Hello guys. Using Maple, I summed the first 200,000 terms and graphed 
this partial sum from x = to x = 50, 000. / also graphed, in yellow, Cy/x, 
where C = J °° u~ 2 smu 2 du (accurate to 8 digits). Of course, in principal, 
the computer could be making all sorts of round off errors, but I thought 
I'd take a look anyway. With all the averaging out that's going on, maybe 
the roundoff errors aren't significant anyway. I've attached the picture to 
this email. To me, the picture suggests that C\fx, [with] C from above, is 
the right thing to try to prove. I agree with you, Michael, that a good way 
to proceed might be by breaking up the sum into two pieces, the 2nd of which 
can be approximated by the integral. 

■■) 

~ Jared" 



7 For convenience of the reader I use M^X to display formulas which Jared described in 
his email. I deliberately left the amusing typo (which happens if you work until midnight!) 
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Here is the picture attached to his email: 
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Fig. 8. (Color online) Speck's 200,000-th partial sum of S2(%) together with a 2 y/x. 

Figure 8 is quite remarkable. The agreement of the displayed graph of 52 (t) 
with that of a 2 Vt, for a 2 given in (JSJ), is nothing less than phenomenal. The 
conjecture that for all t > one has — ce 2 \/i+ "small" fluctuations, with 
a 2 given in (jSJ), had to be right! 

Incidentally, note that the graph of S2 (i) shown in Fig. 8 displays some 
intriguing intermittency, as known from some turbulent flows. For most t the 
graph of S^it) is barely distinguishable from that of a 2 Vi, but every once in a 
while an outburst of activity is visible. How fascinating! Alas, while preparing 
for this presentation, when I replotted the graph with much higher resolution 
in Fig. 2 the intriguing intermittency disappeared. Fig. 8 reinforced Jared's 
intuition about the coefficient a 2 , which would soon be vindicated, but it was 
quite misleading as a guide for how to think about the fluctuations! 

At noon the next day, 52 (t) was the topic of the lunch conversation. In 
particular, Mikko Stenlund, at the time postdoc in our mathphys group, fell 
under the spell of the problem. Two hours later he sent me the following email: 

"Hi, Michael. 

For your information, if instead of the series Yl n sin(x/n 2 ) one considers 
the corresponding integral, one gets with the aid of Fresnel integrals that 



the asymptotic form is a/vtx/2 + sin(s/7r 2 
Mikko" 



sin (x) 



The conjectured a 2 integral gives an elementary value for a 2 — how wonderful! 
Note, though, that the indicated asymptotic expansion (replacing x — > t) is 
for sin(u^ 2 t)dt = y/i J Q /i~ 2 sin(/i 2 )<i/z, not the integral at r.h.s.ffTUl). 

Hardly two hours later Jared had completed our proof of the coefficient 
a 2 . A little upgrading, and also the conjecture about the overall shape of 52 (t) 
was proved. Our proof easily generalizes to all p>l, to which I turn next. 
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3.2 The overall trend of S p (t) 

Analogous to the reasoning for when p = 2, with t > 0, we now split the 
summation in the series defining S p (t) at n = [(pt/r) 1 ^^] =: N p (t/r), thus 

N p (t/r) oo 

Sp(t) = £ sm(n- p t) + £ sin(n-Pt). (11) 

n=l n=JV p (Vr)+l 

When r is small enough (again t < n/2 will do when t gets large), then for 
n > \{pt/r) 1 ^ p+l ' >1 \ any two consecutive arguments t/n p and t/(n + l) p of the 
sine functions will come to lie within one quarter period of sine. Moreover, 
with increasing n, for fixed t/r, the consecutive arguments t/n 2 and tj (n+ l) 2 
will be more and more closely spaced. In other words, for fixed sufficiently 
small r, with increasing t the part of the sum of S p {t) with n > N p (t/r) 
will be an increasingly better Riemann sum approximation of the integral 
In (t/r)+i smiv'Ptydv. Thus, and after the variable substitution v~ p t = £, 

r t/(N p (t/r)+iy 

£ sm(n- p t) w t 1/p i / r 1_1/p sin^. (12) 

Since p > 1, the upper limit of integration at r.h.s.flT2|) goes to oo like At 1 ^^ 
when t — )• oo, and the limiting integral can be evaluated by contour integration: 

POO 

I r^sin^^r^-ijsin^). (13) 

J 

Remark: Integral ( JTBl is related by variable substitution to the generalized 
Fresnel integral J °° sm(r] q )dr] = r(l + ^) sin (f-), which converges for |g| > 1. 

We will sharpen "w ..." in ( I12p to "= ...+ a subdominant error bound." 
This, a similar estimate comparing r.h.s. (p|) with t^x r.h.s. (JT3J) , and the 
subdominance of the first sum in (ITT]) compared to r.h.s. ( IT2|) . leads to: 

Theorem 1. For all p > 1, and all t 6E, we have 

S p (t) = a p sign(t)\t\ 1 / p + 0(\t\ 1 ^), (14) 
wift a p given by r.h.s. (fi^). 

Proof: By the anti-symmetry of 5 p (t) it suffices to consider t > 0, though 
we need to distinguish t < t p and t > t p for some t p > 0. Recall that 
f(pt/r) 1/(p+1) ] =: JVp(t/r). In all estimates below, C is a generic constant. 

First of all, for t p > sufficiently small, we have = At + 0(t 3 ) for 
t < t p , so obviously |5 p (t) - a p t l / p \ < Ct 1/ip+1) for some C when t < t p . 
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Turning to t > t p , for the first sum at r.h.s.(ITT]) the triangle inequality and 
then | sin£| < 1, summing, and an obvious estimate, yield 



N p (t/r) 

^2 sin(n p t) 

n=l 



< 



For the second sum at r.h.s. ffTTj ) we find (for some v n e [n, n 

sm(v~ p t)dv = 



(15) 



1]) 



sm(n p t) 

n=Np(t/r)+l 



( sin(n-Pt) 

n=N p (t/r)+l ^ 

oo 



N p {t/r)+l 
n+l 



sin(t/ p t)dv 

OO / 

( sm(n~ p t) - sm(u- p t) 



n=N p (t/r)+ 



E 



n=N p (t/r)+l Jt/vl 
oo ft/n p 

E / 

n=N p (t/r)+l Jt/uK 



cos£<i£ 



< 



COS £ I dt; < 



E < 

n=AT p (t/r)+l 
oo 

E H^y ~ ( n +l)p) ~~ 
n=AT p (t/r)+l 

t\(pt/r) 1/{p+l) + l\- p < Ct 1/{p+1 \ 



(16) 
(17) 
(18) 

(19) 

(20) 

(21) 

(22) 
(23) 



In this stringjf] as for the first three equalities: ( f!6l) is manifestly true, whereas 
( I17p holds by the mean value theorem for some v n e [n, n+ 1], and (f!8l) holds 
by the fundamental theorem of calculus; as for the ensuing three inequalities: 
f fT9|) holds by the triangle inequality, f )20|) holds since | cos£| < 1, followed by 
elementary integration, while (1211) is due to the monotonic decrease of v t- >■ z^~ p 
for p > 1, with z/ n e [n, n + 1]; the ensuing equality f l2"2"j) holds because the 
sum at l.h.s.f l2"2"|) is telescoping; lastly, inequality ( }2"5|) is obvious. 
For the integral in (Tl6|) the variable substitution v~ p t = £ yields 



p 



t/(Ar p (t/r)+l)f 

r 1 " 1/p sin^ = t 1 / p 



r 1_1/p sin^^ . (24) 

t/(N p (t/r)+l)P J 



8 I am heeding the advice Michael Fisher gave me (after reading [KieOSj ) on Dec. 14, 2007: 
11 1 have now had a chance to delve further into your write-up. Eventually, I found out why 
you say "mean field". The answer is three totally unnumbered equations: That represents 
very bad practice! [...] Please do number all crucial equations in your future papers!" 
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Using one last time the triangle inequality and | sin£| < 1, we find (for t > 1): 



t i/Pl 
v 



-1-1/p 



sin£d£ 



< [(pt/r) 1/(p+1) + 1] < (25) 



lt/(N p (t/r)+l)P 

The entirely elementary proof of Theorem 1 is complete. □ 

Thm.l is illustrated below by three graphs of S p (t) together with their 
trends a p t 1 / p , forp = 3/2,p = 2, and p = y/7. The t interval is always [0, 600]. 
We begin with p = 2 and p = y/7, shown together in Fig. 9. 



30 



20 



io 




500 



600 



Fig. 9. The 5,000-th partial sums of S p (t) for p = 2 and p = ^/7, together with 
their trend functions and T(l — -^=) sin (^y)* 1 ^^: respectively. 

The case p = 3/2, shown in Fig. 10, is interesting in its own right: 

160 




5 00 



600 



Fig.10. The 300,000-th partial sum of 5 3 / 2 (t) together with T(i) sin (f )t 2/3 . 
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Remarkably, a "staircase" structure is clearly visible in the graph of S3/2 (t) over 
the t-interval [0,200], after which it gets more "noisy," yet for 500 <t< 550 
another plateau shows. Doesn't this call for a number-theoretical explanation? 

Since for a moderately small p value like 1.5 a very large number of terms in 
the partial sum of S3/2 (t) was required to achieve a decently converged result, 
I didn't try to push p close to 1; except, a mildly smaller, irrational p = y/2 
was chosen for Fig.6, with a similar expenditure in mode numbers. 

In all three cases shown, the trend function a t 1 ^ truly traces the visi- 
ble trend of S p (t). The erratic fluctuations about the trend are more slowly 
growing in amplitude than the trend. Our Thm.l says that they are bounded 
in amplitude by O^ 1 /^ 1 )). To get an idea of how accurate this bound is, I 
resorted to Maple to plot S p (t) - a p t 1/p =: AS p (t) together with i/yVfr+i) for 
p = 2 and p = y/7, with empirically near-optimized (3 p , see Figs. 11 and 12: 




o 100 200 300 400 500 

t 




O IOO 200 300 400 500 

t 

Fig.12. The 20,000-th partial sum of S^ 7 (t) - T(l - ) sin (^77)t 1/VY with ±0.77t 1 ^ 1 +^ 7 \ 
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Figs. 11 and 12 reveal that the functions ±/3 p t 1 ^ p+1 - ) are accurately bounding 
the growth of the largest fluctuation amplitudes over the shown t interval with 
empirically optimized (3 P ; I have not tried analytically to optimize /3 P . Of 
course, a larger sample of p values would allow a more reliable conclusion; 
however, since much higher precision was needed for these figures, I plotted 
only the cases p = 2 and p = y/7. In section 4.2 we will use a change of 
variables which allows us to plot As p for larger t values, see Figs. 14 & 15. 

After this three-day flurry of activity the inquiry into S p (t) stopped almost 
as abruptly as it had started, because all participants had to return to their 
own important businesses. However, pre-conditioned by my upbringing in 
statistical mechanics, I resolved to resume the inquiry into the fluctuations of 
S p (t) — a p t l f p whenever the opportunity would arise. 



4 Statistics of the fluctuations of S p (t) - a/'P 
4.1 Kac's central limit theorem 

Years later, in March 2011, while listening to Felix Izrailev's interesting pre- 
sentation about quantum thermalization at Michael Kastner's STlaS workshop 
in Stellenbosch, South Africa, I noticed that he refered to work by Mark Kac 
on the central limit theorem for certain trigonometric series. I immediately 
wondered whether this was the information I had been waiting for to hear! 

The relevant original publications are |Kac38] and |Kac 43j , which together 
with some other works by Kac were expanded into his book [Kac59] . According 
to the charming Kac memoir by Henry McKean (cf. p. 219 in [McK90] iR. in 
[Kac38j and [Kac43j the following is proved: 

Theorem 2. Let the set of frequencies {u n } n ^ be linearly independent over 
Q (i.e., for any N G N, the only solution to Yln=i z n u n = with all z„GZ is 
Z\ — ■ ■ ■ — Zn — 0)- Let "meas" denote Lebesgue measure on K. Then 

lim Km ^measjt e [0, T] : a < £ shiest <b} = ^J e^dy. (26) 

Several remarks are in order: in |Kac43j . and in more detail again in 
[Kac59] . Kac himself derives a similar formula in which cosines replace the 
sine functions, and with the t-average taken over the interval [-T, T] rather 
than [0, T\. For Kac's cosine theorem it is obvious (since cosine is an even 
function) that the average over [— T, T] equals the one over [0, T]; however, 
this is not true for a sum of odd sine functions. One has to go through Kac's 
cosines proof to see that, after replacing [— T, T] with [0, T] averages, one can 
also work with sine replacing cosine, indeed. 



3 Note a typo 1/2, instead of l/\/2, in the pertinent formula on p. 219 in |McK90] . 
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At first glance Kac's Theorem looks just like "what the doctor ordered" for 
Spit). Unfortunately, what it says about u S p (t) n (for suitable p) is not about 
S p (t) — instead, it is about the infinite time averages of the family of partial 
sums of S p (t). Of course, S p (t) is defined as the limit iV — > oo of the sequence of 
its iV-th partial sums, but this limit does not commute with the limit T — > oo 
of time-averages over intervals of length T; only for a fixed partial sum, t- 
averaging over [0,T] and summation do commute. Kac's theorem demands 
that for any iV-th partial sum of S p {t) one first performs the limit T — > oo for 
the average amount of time this partial sum eventually spends in the interval 
[ay/N/2,by/N/2], then lets iV -»■ oo (cf. [Kac43j . |Kac59j ). 

Kac's theorem implements Steinhaus' notion of "statistical independence of 
functions:" the average amounts of time which individual sine functions with 
incommensurate frequencies spend in any infinitesimal interval within [—1,1] 
are eventually i.i.d. random variables with mean zero and standard deviation 
1 / y/2 — thus the central limit theorem type appearance of his theorem. Recall 
that the limit T — > oo of the unrestricted t-average over [0, T] of each sine 
function vanishes whereas the t-average of its square converges to 1/2. 

By contrast, we need to take a time average after the infinite summation 
over all sine functions has been carried out and the trend function subtracted. 
This makes it plain that Theorem 2 above is not applicable to our problem! 

Now, all this does not mean that the fluctuations of S p (t) — a p t l l p are not 
normal — they may well be ( "Not so!" is the comment of one of the referees - 
see below). At the time of my SMM 106 presentation, under the spell of Kac's 
central limit theorem, I indeed conjectured that, after "suitable p-dependent 
rescaling," a normal law should hold for the fluctuations of S p (t) — apt 1 ' 13 , at 
least for irrational p (NB: As explained in the introduction, another referee 
noted that p ^ Q won't be sufficient to guarantee the rational linear inde- 
pendence of the involved frequencies {n~ p } ne N-) More precisely, a careful in- 
spection of Kac's proof, which is based on Levy's rigorous version of Markov's 
method of characteristic functions (i.e., Fourier transforms of probability mea- 
sures), reveals that Markov's method should also determine the distribution of 
the fluctuating values of S p (t) — a p t 1//p , yet it also is clear that from some point 
on Kac's arguments will have to be modified. More to the point, even though 
our theoretical error bounds are too rough to show it, empirically the second 
term at r.h.s. fllip seems to be a very accurate Riemann sum approximation 
for the trend function when t/r becomes moderately large, which means that 
As p {t) should be well approximated by the first term at r.h.s. ffTTj ). which is a 
finite sum at each t, containing not more than N p (T/t) terms in the t-averages 
over [0, T\. Since N p (T/r) x CT l ^ p+l \ the longest wavelength in the partial 
sum of sines, which is averaged over [0, T], grows basically oc T p ^ p+1 \ i.e. sub- 
linear in T so that, as T grows large, even the sine functions with the longest 
wavelengths in the partial sum are averaged over many cycles, infinitely many 
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in the limit T — > oo. The upshot is the following conjecture (extending Kac's 
"central limit theorem"): 

Conjecture 1. Suppose p is chosen so that the frequencies {n~ p } ne ^ are ra- 
tionally linear independent. Then 

r , , r b l a 

lmr^measjt G [0,T\: a < ^j^t) £ sin(rr p t) < bj = ^ j e~v y dy. (27) 

Note that Conjecture 1 is weaker than conjecturing that the fluctuations of 
S p (t) — a p t l l p = As p (t) itself are eventually normal. In any event, Conjecture 1 
is interesting in its own right and, if true, may serve as an important stepping 
stone on the way to characterizing the fluctuations of As p (t). I hope to settle 
this issue at some point, or that someone else will feel inspired to do so. 



4.1.1 Experts' opinions on the fluctuations of S p (t) — a p t l l p 

All three referees confirmed my hunch that the questions raised above, and 
many more, can be answered with the techniques of analytical number theory, 
see the books |Vin04j . |GrKo9lj . and the surveys |Mon94j . |IwKo 04j . In par- 
ticular, one of the referees noted: "Assuming the Riemann hypothesis, which 
is quite normal in such type [of] studies, I was able to show that for any e > 0, 

Sp(t) = a p t 1/p + t 1/2(p+1)+£ ^, £ (lnt), (28) 

where g PtE is an L 2 -function, so that \g Pj£ (u)\ 2 du < oo. I think that this 
gives a much stronger control of the error term than [ (1T4"1) ]." Another referee 
made a similar observation about S p (t) as expressed in fl28|) . namely that by 
standard techniques of Fourier analysis and analytic number theory as ex- 
plained, e.g., in [GrKo91j . one would be able to prove that 

Sp (t) = a/'? + t^^Epit); t > 1, (29) 

where E p {t) = 0(t°^) (for t > 1). This referee further noted that "at least in 
some cases (like p = 2) it seems possible to use the methods from the paper 
[HcBr92j to understand the distribution of E p (t). In any case, I think it is easy 
to show that for any p the distribution [of E p (t)] is not going to be normal 
... which would show that this conjecture [about the fluctuations of As p (t)] 
made by the author is wrong." (The referee also outlined the exponential sum 
obtained for E p (t) with the methods in |GrKo9l] . and why the distribution of 
E p (t) should be non- normal. Unfortunately, since I am no expert in analytical 
number theory, I refrain from making an (inevitably amateurish) attempt to 
explain these arguments here.) The third referee similarly noted that "the 
Gaussian behavior is violated in other examples in number theory under a 
similar philosophy |BCDL93j (see also |Ble92] V 
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I am grateful to stand rectified about my speculations about the normality 
of the fluctuations of As p (t). I am also grateful for the observation that p £ Q 
is insufficient to guarantee rational linear independency of the set of frequencies 
{n~ p }. As far as I can see, though, my (so revised) Conjecture 1 about the 
normality of fluctuations of the [t 1 /( p+1 ']-th partial sum of S p (t) is still viable. 

4.2 Late-t asymptotics of S p (t): Riemann's £ function 

Sometimes wondrous things happen. After reading the announcements of the 
SMM 106 conference talks, Norm Frankel and Steve Miller requested the pdf 
file of my upcoming talk. Since, as usual, I finished the preparations for my 
talk barely in time, they had to wait until then. Soon after, they got back to 
me with exciting emails, the essence of each of which I pool together. Here 
first is Norm Frankel (whose inimitable style I like to preserve): 

"Dear Michael, 
MILLE GRAZIE! 

I've been looking forward to receiving this and will read and study it with 
relish [and mustard - HI]. I'll write back when I have > ln(2) to say. . . . 
Using a Mellin transform, the asymptotics comes out in one line. I find 
S2{x) ~ (1/4)a/27t/x — 7r/ 4+ intricate terms. Similarly S p (x) can be readily 
exhibited. I seem to be differing by a factor of 1/2. . . . 
The results I sent were for 5_ 2 (x), not Si{x). OF COURSE yours is a 
MUCH trickier sum — back 'gain soon. . . . 

I've just looked again at your sum. I think what I started out to do is 
correct in concert with analytical continuation: the Mellin transform of 
sm(ri p x) = sm(irs/2)r[s]((ps)x~ s . Inverting readily gives the large x 
asymptotics — even for p = —2, your series — correction terms follow. 
It may be incorrect; if so, mia culps. . . . 
THE HAPPIEST OF HOLIDAYS FOR YOU AND YOURS! 
Warmest Regards, Norm" 
And here is Steve Miller: 
"Dear Michael, 

thank you very much for your e-mail. You really made excellent powerpoint 
slides - 1 am sure the audience appreciated them. They were very clear and 
entertaining, [...ffl I had not seen anything with n~ 2 like you have. 
I have a few comments which I hope will be helpful to you. You asked about 
the connection with the Riemann ( function at the end of your talk. If you 
(formally, at least) take the Mellin transform of the function f:y*-^-sm(x/y 2 ) 
which is Mf(s) = — ^x s ^ 2 T(— s/2) sin(7rs/4), then your sum is a contour 
integral of Mf(s) times the Riemann ( function. The Riemann ( function 
has only one pole, at s = 1, and the residue of Mf(s)((s) at s = 1 is your 

10 Steve drew my attention to the intriguing papers |ChUb07l IMiSch04] on sin(ro 2 :c) series. 
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main term a/vtx/2. The only other poles of M f(s)((s) are at values of s 
of the form 2 + An, n > 0/ these give correction terms in the asymptotic 
expansion of Greenfield's infinite sum as x — > 0. So the full asymptotic 
expansion should come from this. 

I also want to note that the functional equation of the Riemann ( function 
(again, completely formally - this is essentially Poisson summation here) 
gives the identity that your Greenfield's sum is the sum of g(n,x) over 
n > 0, where the Mellin transform of g(z,x) in z is 

-(2tt) 1 - s x( 1 - s )/ 2 cos(tt( S + l)/4)r(s - l)C(s)/r((l - s)/2) 
This function (without the ((s) factor) is a sum of hypergeometric functions 
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Fig. 13. The sum of hypergeometric functions Steve included as image. 

Mellin inversion shows Greenfield's sum is also a sum of this function. 
That dual point of view may reveal some other properties of the sums you 
considered in your slides. [I think I botched a calculation here, but I hope 
the idea and strategy was clear]. 
Best holiday wishes, 
Steve" 

I was thrilled — and immediately turned contemplative: While Norm is 
located half around the globe away from my office, Steve Miller's office isn't 
much further away from mine than Steve Greenfield's! Here I was, having the 
problem in the back of my mind all these years — without ever mentioning 
it to Steve (M.)? What if Steve G. would have sent his question to Steve M. 
instead? It reminded me of "Missed Opportunities" |Dys72| , the beautiful 
Gibbs lecture by Norm Frankel's longtime friend Freeman Dyson. 

Enlightened by their comments I decided to compute the late time asymp- 
totics of Spit) beyond the leading order term. However, things aren't quite as 
straightforward as they seem! 

First, let me flesh out what Norm and Steve wrote in their emails. Recall 
that the Mellin transform of a continuous function / : M + — > R is 

/■oo 

(Mf)(s) := / y s - l f(y)dy (30) 

wherever r.h.s.( l30i) is well-defined. Supposing \f(y)\ = 0{y~ a ) for y J, and 
\f(y)\ = 0{y~ b ) for y •f oo, with a < b sharp, then (M/)(s) is analytic in its 
fundamental strip a < < b, where it tends to zero as \$s(s)\ —> oo (by the 
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Riemann-Lebesgue lemma). Writing (M/)(s) =: f(s), its inverse transform is 
given by the straight contour integral (in the improper Riemann sense) 

rc+ioo 

f(y) = h \ y- s f(s)ds any c £ (a, b). (31) 

J c—ioo 

Turning to S p (t), for t > we define r := t~ l l p and write 5 p (t) = %p(r), i.e. 

= E„ eN sin ((nr)-P); p > 1. (32) 

Since |^>(r)| = 0(r _1 ) for r | and |^>(r)| = 0(r~ p ) forr f oo, for 1 < ?R,s < p 
we can take the Mellin transform of %,(r) to find, with obvious manipulations, 

POO 

%(s) = / r 3 - 1 sin ((nr)- p )dr (33) 

J0 raeN 
poo 

= E / r^ 1 sin ((nr)~ p )dr (34) 

n6NJ0 

/■oo 

= E£ / y^sin^)^ (35) 

neN JO 

= C(*)|r(- f)sin(- f{). (36) 
which is analytic in 1 < 9fts < p, its fundamental strip. And so, for c £ 

/•e+ioo 

= 2^ / ^C( S )ir(- J) sin(- f j)<fc, (37) 

or, after switching back to t — r~ p , renaming s/p into q, and introducing 

Q,(?) = C(P?)r(-?)sni(-|?), (38) 
we find 

5 p (t) = JL / f^( f )cfc. (39) 

•/ c/p—ioo 

This is how far you can get using only Euler's series for the Riemann ( function. 

Next we use that r.h.s. (l36l) . understood as analytic extension of l.h.s. f l33j) . 
is manifestly meromorphic in C, having simple poles at s = 1 (coming from 
the Riemann ( function) and at s/p = 2n — 1, n £ N (coming from the Euler 
T function); the poles of the T function at s/p = 2(n — 1), n £ N, are ironed 
out by the pertinent zeros of the sine function. Note, though, that the ( pole 
and the T poles are located on different sides of the fundamental strip. 

Therefore, if we now shift the contour in the s plane to the right, beyond 
all T poles, we obtain the Taylor series expansion of S p (t) about t — 0, viz. 

OO 

s P (t) = E (-i) & ^%SS^ 2fe+1 ; ®p > i- (40) 



k=0 
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It is readily checked that the same expansion is obtained directly from ([I]) by 
replacing sin(n~ p t) by its Maclaurin expansion, then exchanging the Maclaurin 
summation with the summation over nGN given in and using the Euler 
series of the Riemann ( function for s > 1. 

If, on the other hand, we shift the contour to the left just a little bit beyond 
the pole at s — 1, say to c p = 1 — e, we pick up the pole's residue and obtain 

p — +ioo 

^(t) = a/'? + JL / " f Q,(0<fc, (41) 

J — — ioo 

v 

with a p given in ( 1T31) . Incidentally, by corollary to Theorem 1, the integral at 
r.h.s. (l4"Tj ) is bounded in magnitude by (3 p \t\ l ^ p+l \ 

So we see how the Mellin transform plus the residue theorem of complex 
analysis reproduces — in one elegant sweep — all the results we could establish 
with more elementary means, save the bound on the integral at r.h.s. fj4Tj ). 
Alas, with (T4"0"j) and ( 14 ip we exhaust the information about S p (t) which one 
can extract from the poles of Qj(<?). Clearly I cannot end on such a note! 

Let's see what we can learn from the fact that the fluctuations As p (t) about 
the trend ct p t 1//p are given by the contour integral at r.h.s.Pdj). Proceeding now 
first formally, we pretend that we can shift the contour in f l4~Tj) to any c p < 1. 
Writing it as {<; = c p /p+iv', v' e M} and then changing the integration variable 
in the contour integral at r.h.s.(l4"Tj) to v', and then to v — v' + ic p /p, gives 

/■Cp/p+ioo 

A 4(*) = h\ * ( 42 ) 

J Cp/p—ioo 

= t c ^(5- 1 Q p (^ + ^-))(lnt) (43) 
= (?-%(i-)){]nt), (44) 



where 



/oo 
e™f(v)dv (45) 
-oo 



denotes the formal inverse (non- unitary) Fourier transform of f{v). To make 
proper sense out of the formal integral (HSJ), Fourier analysis enters next. 

Since the appropriate Fourier variable for As p (t), t > 0, is u = hit, Figs. 11 
and 12 are a bit deceptive now. The graph of u (-> As p (e u ), «eR, shown 
below for p = 2, is a better guide to one's intuition. The u interval corresponds 
roughly to the t interval in Fig. 11, though not quite (note that e 6 ~ 400). 
Fig. 14 reveals that A52(e") has only one zero to the left of u = and vanishes 
x — a 2 e M / 2 when a | — oo. On the other side of u = 0, A52(e n ) develops 
oscillations with wavelengths which become exponentially small as u gets large 
while their amplitudes grow with u, though bounded by the theoretical bounds 
±/3 2 e M//3 , with empirically optimized /?2 = 20/29 (cf. Fig. 11). 
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Fig.14. The graph ofuK A5 2 (e"), u £ (-4, 8), together with ±|§ exp(u/3). 



What Fig.14 only hints at is that beyond u = 5.5 the growth of the ampli- 
tudes departs more and more from our theoretical bound. This is illustrated in 
Fig. 15, which is the continuation of Fig.14 to the right — rescaled, of course, 
to fit on this page. 
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Fig.15. The graph of u M- AS2(e u ), u e (8, 20), together with ±|§ exp(u/3). 



Fig.15 leaves no doubt that the bounds ±/3 2 exp(w/3) become lousy for large 
u; in fact, something like ±72exp(w/5) traces the fluctuation amplitudes much 
better for u G (8,20), but for even larger u (not shown) also this bound will 
outgrow the fluctuations. (NB: According to the referees, see 4.1.1, the correct 
bound should be ±k 2 exp(w/ (6 + o(l))).) Nevertheless, as Fig.15 indicates, the 
fluctuations continue to grow forever. To prove their unbounded growth is not 
so easy, but at least it is readily shown, for all p > 1, that As p (e u ) does not 
approach when u — > oo — for suppose it would, then also As p (t) — > when 
i | oo, and so then does its t derivative (because t^ p ^ p I when t | oo, 
and Spit) contains a smallest wavelength); but the t derivative of S p (t) is a 
manifestly quasi-periodic function of t: a contradiction — end of proof. 
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The upshot of this discussion is that u >->■ As p (e u ) is a tempered distribu- 
tion. Therefore, its Fourier transform v i-> (^[Aj^ o exp])(t> ) = Qp(iv) is to be 
understood in the sense of tempered distributions as well. 

In this vein, let & denote the Schwartz space of complex £°° functions on 
R which together with all their derivatives decay to zero at infinity faster than 
any power. If ip G &, then its Fourier transform Jip G &, too, where 



= / e- luv i)(u)du. (46) 
Jr 

The Fourier transform of a tempered distribution g G & is then defined by 
(3 r ^)(v)(3 r "V)(^)^ = / g(u)i;(u)du V ^ G ©, (47) 



where I hope to be forgiven for using the merely formal integral notation rather 
than a proper dual pairing notation, cf. [ReSi75] . 

As to the real function u h-> Aj> p (e u ), for our purposes it suffices to inspect 
its properties when integrated against the members of the family of shifted, 
scaled Hermite functions {^ n (/c(u — w)) G & : w G R, k G M + }^L , with 

V„(w) = (2^! v ^)- 1 / 2 e-" 2 / 2 ^ n ( M ), (48) 

where 

= (-l) n e u2 £re- 2 (49) 

is the n-th Hermite polynomial, with Hq = 1. In particular, to determine the 
late t asymptotics of As p (t), we now define quantities of the form 

A p n (w;n):= [ As p {e u )^ n {n{u - w))du (50) 



and evaluate their asymptotics as w — > oo with the help of (|4j 

Recalling that the Fourier transform (^[ASp o exp])(t>) = Q p {iv) 1 we have 



A>(w;k) = / ^(^[io^. -w)])(v)dv. (51) 

The integrand can be recast into a more convenient format by noting that 

(^„o(«(. -w))){v) = j™\(3-^ n )(^ (52) 
and by recalling that the Hermite functions are £ 2 eigenfunctions for 3 r , vizi"] 

(^V»)(«) = ^»(«)- (53) 



11 The factor is a consequence of working with the non- unitary version of 3 ". 
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Note that v h-> ip n {v) is even for even n and odd for odd n. Furthermore, 
consulting |Edw74] (or |Tit86J . Jlvi03j ). one sees that it follows directly from 
the explicit formula (1381) thatP*! v i— y 3?(Qp(if)) is even and v (->■ S(Qp(if)) is 
odd, shown for p = 2 in Fig. 16: 




-lO -5 O 5 lO 



Fig. 16. Real (solid) and imaginary (dashed) parts of v t— > 02(^)1 f S K. 

Therefore, the real part of the function v H- e mw Qp(iv) is even, while its 
imaginary part is odd. In summary, we can conclude that (l5"Tj) simplifies to 

POO 

A p n {w-K) = {-ifiJm ^(D^iv^dv (n even), (54) 

v Jo 

resp. 

POO 

Aliw-n) = (-l^JlQ ty n (l)Qe{iv)e ivw dv (n odd). (55) 

v Jo 

The pertinent real or imaginary part of the integrand which features in the 
integrals at r.h.s.( ]5"4"|) . (l5"5"|) is shown for n = and n = 1 in Figs. 17 & 18 below, 
respectively, in each case with p = 2, k = 1, and u> = 60. 
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Fig.17. The graph of ?R(e ivw Q2(iv)) tp (v) for w = 60. 



12 The Fourier transform v H> Qp(iv) of the real function w i— > As p (e u ) cannot be purely 
real or purely imaginary, for As p (e u ) is neither even nor odd, see Fig. 14. 
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Fig.18. The graph of ^{e ivw Qz(iv))ili\{v) for w = 60. 

Since for each choice of n, p, k the integrand of the integrals at r.h.s. fp>4"]) . (l5"5"]) 
becomes a highly oscillatory Schwartz function, by the Riemann-Lebesgue 
lemma the A^(w; k) vanish in the limit w — > oo, and we are interested in how 
they vanish asymptotically when w — > oo. Since the integrand is a complex 
analytic function in a strip neighborhood of R, their w — > oo asymptotic 
behavior is easily found from the full asymptotic expansion 

poo l 

/ iu~)^y™<iv = e i[^(a^H](o) e ^m^-m 

Jo eeN j=o 

(Here, [•••] does not mean "integer part.") Inserting ( 15T)|) at r.h.s. (l51]) . (|B3]) 
yields the aysmptotic w — > oo expansion of A^{w; k). For most practical 
matters we will only need the pertinent lowest order nonvanishing term in 
( )56|) . as in the example below. 

4.2.1 Application: Gaussian averages of u i— > As p (e u ) 

I briefly register the bonus of the asymptotics obtained in the previous sub- 
section: since ipo(u) = ix~ x ^e~ u / 2 , it follows that Aq(w; 1) is proportional to 
a standard Gaussian average of As p (e u ) centered at w. For instance, 

^;l) = ^(|T + ln(27r))i + 0(^), (57) 

where 7 = 0.57721... is Euler's constant. Multiplication of ( 1571) by 7r 1//4 /v / 2~7r 
yields the standard Gaussian average of Aj> 2 (e") centered at w. 

The road is now paved to analyze the fluctuations of 11 4 As p {e u ) as 
w 00. If one thinks of the w-centered Gaussian u- average as an analogue of 
Kac's uniform t-average, one can ask how much "u-time" the function As p (e u ) 
spends on average, centered at w, in some value interval (a,b). Markov's 
method can be applied, and one might even be able to compute the answer 
explicitly using the formalism of the previous subsection, whether the answer 
is "normal" or not. So much for the Mellin transform. 
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5 Epilogue 



After my presentation, which ended with sect. 4.1, a young participant at the 
meeting came to talk to me. He was incredulous, and something close to the 
following conversation ensued 1^1 
"Prof. Kiessling: 

Why are you doing this? What does this have to do with physics?" 

"Well, I am doing this because it's fun! And didn't I myself at the end of 
the talk raise the question whether there are connections to physics? Also, 
recall that I noted that the Schrodinger hydrogen spectrum has eigenvalues 
— n~ 2 , so n~ 2 frequencies do occur in physics, you get cos(n _2 t)-M sin(n~ 2 t) 
as t-dependent factors in any expansion using the eigenwave functions. Of 
course, the series S2{t) itself may not occur." 

"Exactly, how can you hope that somehow this will be useful to physics? 
This is crazy!" 

Truth be told, I am not sure he really said "crazy," but he certainly gave 
me the impression that the thought crossed his mind. And on this anecdote 
I close by quoting Jerry Percus from an interview he gave a few years ago at 
the Courant Institute [Bal08j : 

"What you want is to be a little bit crazy. You want to think of things 
that sound like nonsense to start with and then when you get deeper, 
they're not nonsense at all." 

**# 

Acknowledgement: I thank Steve Greenfield for his original question which 
triggered this line of inquiry; I also wish him a happy retirement. Next I thank 
Jared Speck and Mikko Stenlund for their participation in nailing down a 2 - 
Thanks go to Norm Frankel and Steve Miller for their enlightening explana- 
tions of the relationship between S p (t) and £(s) after my SMM 106 talk; for 
further enlightenment see also [ NHFG921 IKFGT95] . Of course, I thank all 
of these, and also Michael Fisher, for allowing me to quote from their emails 
they had sent me. I also thank Joel Lebowitz for the honorable invitation to 
contribute to SMM 106 and to this special issue, and the three anonymous 
referees for their generously offered expertise in analytical number theory (in- 
cluding 11 references, amongst them also [HaLi36] IFle50] ) which allowed me 
to revise the perspective on S p (t) presented in this paper; it also allowed me 
to correct my Conjecture 1. But most of all, I thank all three honorees for all 
they taught (all of) us at so many a statistical mechanics meeting, and Jerry 
Percus in particular — for being a wonderful mentor and friend! 

13 Many of the historical developments and anecdotes I wrote up for this talk have been 
faithfully reconstructed from email records, but some of it only from my memory, and so it 
is appropriate to quote another of Jerry's favorites: "Don't mistake me for the facts!" 
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Appendix 

Here we fulfill Greenfield's request and produce lower and upper bounds on 
(i) suitable for an undergraduate workshop. 
First of all, pick iV G N U {0} and split r.h.s.([T|) with p = 2 into two parts, 

N oo 

S 2 (t) = £ sin(rt) + £ sin(n- 2 t), (58) 

n=l n =N+l 

with En=i = 0- Now sin(£) > —1 estimates the first sum from below by — N, 
and sin(£) > £ — £ 3 /6 for £ > is used to bound the second sum by 

oo / oo \ / oo \ 

£ sin(n- 2 t) > ( £ r2- 2 )t-|( £ n- 6 U 3 . (59) 
The converging p-series are easy to evaluate when A" is not too big, using 

oo N 

E n- 2 = ^-5>" 2 , (60) 

n=N+l n=l 

oo N 

E ^- 6 = ^-E^- 6 - (6i) 

n=iV+l n=l 

And so, with the abbreviations r.h.s. fl60]) = A 2 (N) and |x r.h.s. flBT]) = B 2 (N), 
S 2 (t) > -N + A 2 (N)t - B 2 (N)t 3 ViVeNU{0}, (62) 
and thus 

S 2 (t) > max{-N + A 2 (N)t-B 2 (N)t 3 : < N < A*} V A* G NU{0}. (63) 

It suffices to pick A* = 7 to obtain a ready-to-plot small family of cubic 
parabolas, the pointwise maximum of which is a positive lower bound to S 2 {t) 
which "tilts upward" over the whole interval < t < 120. 

For the sake of completeness of the discussion I should note that with equal 
ease one can also produce a complementary upper bound to S 2 (t). Namely, 
instead of sin(£) > — 1 one now uses sin(£) < 1 to estimate the first sum in 
(155|) from above by N, next one uses that sin(£) < £ — £ 3 /6 + £ 5 /120 for £ > 
to estimate the second sum in (I58p from above by 

oo / oo \ / oo \ /oo \ 

E sm(n-H)<( £ n- 2 k-±( E ^V + lsol £ ^V- (64) 

n=JV+l V n=JV+l 7 ^n=7V+l 7 V n=N+l J 

The first two converging p-series are just the same as before, the third one is 
new, but which is equally easy to evaluate when A" is not too big, using 

E ^ = ^-E-- 10 - (65) 

n=N+l n=l 
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And so, with the new abbreviation r.h.s. (l65] ) = C 2 (N), we have 

S 2 (t) < -N + A 2 (N)t-B 2 (N)t 3 + C 2 (N)t 5 ViVeNU {0}, (66) 

and thus, V N* G NU {0}, 

S 2 (t) < min{-iV + A 2 (N)t - B 2 {N)t 3 + C 2 {N)t 5 : < N < N,}. (67) 

Once again picking N* = 7 we now obtain a ready-to-plot small family of quin- 
tic polynomials, the pointwise minimum of which is a positive upper bound to 
S 2 {t) which also "tilts upward" over the whole interval < t < 120. Together 
with the family of lower bounds from above, this does produce an upward tilted 
corridor in which the graph of S 2 {t) must lie. This is illustrated in Fig. 19. 
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Fig.19. The 2,000-th partial sum of S 2 {t) together with the bounds J62]) and $66$ for < N < 7. 
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